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m EXPLORATORY INVESTIGATION OF SOME TYPES OF AEROELASTIC 

INSTABILITY OF OPEN AND CLOSED BODIES OF REVOLUTION 

MOUNTED ON SLENDER STRUTS 

By S. A. Clevenson^ E. Widmayer, Jr., 
and Franklin W. Diederich 


SUMMARY 


Aeroelastlc instability phenomena of isolated open and closed rigid 
bodies of revolution free to move under elastic restraint have been 
investigated experimentally at low speeds by means of models suspended 
at zero angles of attack and yaw on slender flexible struts from a wind- 
tunnel celling. Three types of instability were observed - flutter simi- 
lar to classical bending -tors ion flutter, divergence, and an uncoupled 
oscillatory instability which consists in nonviolent continuous or inter- 
mittent small-amplitude oscillations Involving only angular deformations. 
The speeds at which this oscillatory instability starts were found to be 
as low as about one-third of the speed at flutter or divergence and to 
depend on the shape of the body, particularly that of the afterbody, and 
on the relative location of the elastic axis . 

An attempt has been made to calculate the airspeeds and, in the 
case of the oscillatory phenomena, the frequencies at which these insta- 
bilities occur by using slender-body theory for the aerodynamic forces 
on the bodies and neglecting the aerodynamic forces on the struts. How- 
ever, the agreement between the speeds and frequencies calculated in this 
manner and those actually observed has been found to be generally unsat- 
isfactory; with the exception of the frequencies of the uncoupled oscil- 
lations which could be predicted with fair accuracy. The nature of the 
observed phenomena and of the forces on bodies of revolution suggests 
that a significant improvement in the accuracy of analytical predictions 
of these aeroelastlc instabilities can be had only by taking into account 
the effects of boundary-layer separation on the aerodynamic forces. 
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INTRODUCTION 


Flutter, divergence, and similar aeroelastic instability problems of 
wings and tail surfaces have been recognized for a long time. On the 
other hand, the related problem of aeroelastic instability of bodies of 
revolution {generally hereinafter referred to simply as "bodies") has 
become of interest only recently, primarily because only recently have 
external stores and fuel tanks in the shape of bodies of revolution been 
carried on high-speed airplanes, and only at high speeds do the aero- 
dynamic forces exerted on bodies at low angles of attack become suffi- 
ciently large to give rise to aeroelastic problems. 

There sne several differences in the aeroelastic instability phe- 
nomena of wings and of bodies, that is, in the nature of the motions, in 
the nature of the aerodynamic forces involved, and in the nature of the 
resulting phenomena. 

The aeroelastic phenomena of wings essentially Involve deformations 
of the wings themselves, whereas those of bodies are very unlikely to 
involve significant deformations of the bodies and are based, instead, 
on the deformation of the members supporting the body. For instance, a 
fuel tank carried on two struts, one behind the other, ixnder a wing, or 
a ram jet carried similarly on supports above the fuselage can move later- 
ally as a result of the sidewise deflections of both struts in the same 
direction, and they can be yawed by a deflection of the front strut to 
one side and of the rear strut to the other. In these two degrees of 
freedom, classical flutter may occur under the proper circumstances; 
mder other circumstances and involving only the yawing degree of free- 
dom, classical divergence may occur. 

The aerodynamic forces on wings at small angles of attack or imder- 
going oscillations of small amplitude about zero angle of attack can 
generally be calculated with sufficient accuracy by potential-flow theory; 
they are linear functions of the angle of attack or the amplitude, respec- 
tively, and are not influenced in an essential way by the boundary layer. 
(Exceptions to this statement are the forces causing such nonclassical 
instability phenomena as stall flutter, aileron buzz, and wing biiffeting.) 
The aerodynamic forces on bodies of revolution, however, are often essen- 
tially determined by the effects of viscous flow. For instance, the lift 
which is known to exist on bodies at an angle of attack in steady flow is 
due entirely to these effects, because potential-flow theory predicts zero 
lift for this case. This lift is often an intrinsically nonlinear function 
of the angle of attack. (See ref. 1, for instance.) Consequently, the 
aeroelastic instability phenomena of bodies are more likely to be of a 
nonclassical type related to stall flutter and similar phenomena than are 
the aeroelastic instability phenomena of wings. 
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Two experimental flutter investigations have "been made of bodies 
mounted on wing tips (refs. 2 and but no aeroelastic -instability- 
studies appear to have been made previously of essentially isolated 
bodies, that is, bodies mounted at some distance from a lifting surface 
on flexible struts which contribute no aerodynamic forces. Some instances 
where this problem arises in practice are external stores or tanks car- 
ried on struts \mder the wing and ram Jets carried on supports on top of 
the fuselage. Also, an analysis of the aeroelastic instability of an 
isolated body may serve to shed some li^t on the much more complicated 
problem of aeroelastic instability of a body mounted on a wing tip. 

An investigation has therefore been conducted in order to gain some 
insight into the nature of the instability phenomena of such isolated 
bodies of revolution. A streamlined body, an open tube, and several 
bodies consisting of the tube with -variously shaped end pieces were sus- 
pended from the ceiling of a wind tunnel on struts of several stiffnesses. 
The closed bodies were intended to simulate external stores or fuel tahks| 
the open tube, an unfired ram jet. In one series of tests the tube was 
also mo-unted on two struts covered by a fairing. The nature of the vari- 
ous types of aeroelastic instabili-ty that occtirred under varioxis condi- 
tions was observed, as were the airspeeds at which they occurred and the 
frequencies of any oscillations present. All tests were conducted at low 
speeds (Mach nximbers less than O. 5 ) and -with a range of Reynolds number 

(based on body length) from I .5 X 10^ to 7*1 X 10^. 

In an attempt to analyze some of these results, the speed and oscil- 
latory frequency at which -various types of aeroelas-tic instability may 
occur have been calculated by using slender-body theory for the aero- 
dynamic forces, with certain additional assumptions in the case of the 
open tube. The deri-vation of these forces is presented in the appendix 
of this paper, and the calculations are described therein. The calculated 
and observed results axe compared and discussed. 


SYMBOLS 


a 


f- 


h 


ratio of distance of elastic axis of supporting strut 
system rearward of midpoint of body to one -half length 
of body; in case of bodies consisting of tube with various 


end pieces, the midpoint is that of tube. 


2si 


frequency of oscillatory instability, cps 

natiural lateral bending frequency of body on struts measured 
in still air, cps 
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natural frequency of yawing oscillations of body on struts 
measured in still air, cps 

structural damping coefficient (see, for e xa.Tnp l e, ref. 4) 
lateral translation at elastic axis, ft 
amplitude of lateral-translation oscillation, ft 

moment of inertia of bod.y about elastic axis of configu- 
ration, slug-ft2 

effective lateral bend.ing spring constant of supports with 
body mounted, Ib/ft 

effective yawing spring constant, ft-lb/radian 
reduced- frequency parameter, Lcuy^V 

aerodynamic force per unit length along body, Ib/ft 

length of body (length of tube, in case of bodies con- 
sisting of tube with end pieces), ft 

aerodynamic moment about elastic axis, ft-lb 

mass of body, slugs 

aerodynamic (lateral) force, lb 

dynamic pressure (at onset of instability, unless specified 
otherwise), Ib/sq ft 

dimensionless dynamic -pressure parameter, 2qVb/Kix 

radius of body of revolution, ft 
dimensionless cross-sectional area of body, 

coordinate along length of body, measured rearward from 
nose, ft 

distance from nose of body to elastic axis of support 
system, ft 
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S2 


T 




a 


a 


o 


I 

p 


a 


03 


distance from nose of body to center of gravity of body^ 
ft 

volume of body, cu ft 

speed (at onset of instability, unless otherwise speci- 
fied), ft/sec 

ratio of distance of center gravity of body rearward of 
elastic axis of support system to one-half length of 

body, — -i. 

L/2 

yawing angle, radians 

amplitude of yawing oscillations, radians 

dimensionless coordinate, s/L 
density of test medium, slugs/cu ft 

density of air at standard sea-level condition, 0.002578 
slugs/cu ft 

dimensionless distance from nose of body to elastic axis, 
si/L 

angular frequency of oscillatory instability, radians/sec 

DESCRIPTION OF TESTS 
Apparatus 


¥ind tunnel .- The tests were conducted in air at variable pressures 
in the Langley ^.5 -foot flutter research tunnel, which is of the closed- 
(circular) throat single-return type. 

Models . - The two basic models, the airfoil-shaped body and the open 
tube, are shown in figure l(a) . The airfoil -shaped body of revolution 
is that generated by rotating an airfoil about its chord. The ordinates 
of the airfoil are twice those of an NACA 65 -OIO airfoil and are listed 
in table 1. The open tube consisted of aluminum sheet I /16 inch thick 
rolled into a cylinder with a diameter of 6 inches. Various end sections 
were used in conjiinction with the tube to represent closed bodies of 
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revolution with a cylindrical center sectioni these end sections are 
shown in figure l(h). The open tube and the airfoil-shaped body are shown 
mounted on their supports in the tunnel in figures 2(a) and 2(b)^ respec- 
tively. (The scales shown in fig. 2 read in inches.) 

Strut-support systems .- The models were mounted on one of foirr sets 
of supports which consisted of two small-diameter steel rods fixed on one 
end to a mounting plate, and on the other end to a mounting bar to which 
the models were bolted. (See figs. l(a) and 2.) The struts were designed 
to make the lateral-bending frequency of the bodies about one-half or one- 
third of the frequency of their yawing oscillation in still air. The 
natural frequency of the forward and rearward oscillations of the bodies 
was approximately six times their lateral -bending frequency. The strut 
diameters and effective spring constants of the support systems comprised 
by these struts are listed in table 2. 

The effective spring constants were obtained from the frequencies 
f^ and fo; measured in still air and the known masses and moments of 

Inertia by means of the relations 


Kh = (atfh)^ 

and 

Ka = (2rtfQ:)^lQ: 


The values of and Kq, obtained in this way represent spring constants 

in the true sense of the word only when Xq, = 0, because only then are the 

yawing and sidewise-bending oscillations uncoupled (and even then only if 
the additional apparent mass of the still air is neglected) . The values 
given in table 2 are averages of the values obtained with different bodies 
for Xq; = 0 (except for the values listed for struts A, for which fre- 
quency measurements were made for Xj^ = 0.10 and 0.l4). 

For one series of tests with the open tube, the struts B were covered 

with a fairing of aluminum alloy l/32 inch thick, which extended about 

1 inch ahead of the front strut and inch behind the rear strut and was 

2 

attached to the struts along their entire lengths but was not attached to 
the mounting plate or mounting bar . The airfoil obtained in this manner 
was about 5/l6 inch thick at and between the two struts, had a rounded 
leading edge, and a sharp trailing edge. 
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Strain-gage Instrumentation .- The only instr-umentation^ apart from 
the usual Instrumentation required to measure the tunnel speed and den- 
sity^ consisted of electrical-resistance strain gages moixnted on the 
roots of the struts in such, a way as to measure strains due to sidewise 
deflections of the struts. The output signals of the gages on the front 
and the rear struts were amplified separately and fed to two channels of 
a multiple recording oscillographs and also to an oscilloscope for immedi- 
ate visual observation. 


Tests 

General procedure .- The procedure for each test was to increase the 
tunnel speed slowly and at the same time the angle of yaw of the model 
adjusted (by yawing the mounting plate by means of a mechanism outside 
the tunnel) to give zero lift and moment on the body, as indicated by the 
strain-gage outputs. When some ts^e of instability occurred, the strain- 
gage outputs were recorded, the type of instability was noted, the tun- 
nel conditions were observed, and the test was discontinued, except in 
some instances when the instability was not violent and it was desired to 
study it further. 

The model-strut combinations tested in this manner are listed in 
the left half of table 3 } also listed are the model mass, model moment 
of inertia, elastic -axis and center-of -gravity location, as well as the 
measured still-air frequencies of each configinration. The tests are 
divided into several series for the sake of convenience in referring to 
them. 


Tests on the streamlined body .- In series I, the streamlined body 
was moimted on struts A, the most flexible ones, and the tests were con- 
ducted at various air densities. In series II, the same body was momited 
on the somewhat stiffen struts B. Only the nonviolent uncoupled oscil- 
latory instability occurred even at the highest air density used in these 
tests (substantially sea-level density) . The tunnel speed was increased 
successively to several values beyond that at which this instability first 
occurred, the air density being kept substantially constant at the sea- 
level value. Test series II consists of measurements of the frequencies 
of the oscillations at these air speeds. . The streamlined body was also 
mounted on the still stiffen struts C, but no aeroelastic instability of 
any kind occurred; consequently, this experiment is not listed in table 3- 

Tests on the open tube with miscellaneous end sections .- The open 
tube with various end sections was mounted on struts B, C, and D, that 
is, on all but the most flexible struts. Series III consists of the tests 
made with the various configurations at constant air density. Series IV 
oconsists of tests made at various densities by using the tube with hemi- 
spherical sections at both ends mounted an struts C . 
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Tests on the open tube on unfalred struts .- Tests made at various 
air densities on the open tube mounted on the most flexible struts com- 
prise series V, and tests made at constant air densities on the open tube 
with various center-of-gravity positions mounted on the next stiffen 
struts comprise series VI. Ihe center of gravity was varied by attaching 
narrow bands of lead l/l6 inch thick to the inner surface of the tube^ so 
that the mass, the moment of inertia, and, hence, the still-air frequencies 
were changed as well. The open tube was also mounted on struts D, but no 
instability occurred at any speed up to and including the speed at which 
this test was discontinued, namely, 536 feet per second. This test is, 
therefore, not listed in table '!). 

Tests on the open tube on faired struts .- The tests which constitute 
series VII are those made at constant air density on the open tube mounted . 
on struts B with the fairing attached. The location of the center of 
gravity was varied in the same manner as in series VI. 


RESULTS 


Presentation of results .- The results of the tests are presented in 
the right half of table 3 and some of these results are presented in fig- 
ures 3 to 9' 

The speeds listed in table 3 for tests which led to flutter or diver- 
gence are those at which these instabilities first occurred. Similarly, 
the frequencies listed for the tests which led to flutter are those at 
the flutter condition. The air speeds listed for the tests of series IV 
are those observed when the instability first occurred, and the frequen- 
cies are those observed at that speed. Similarly, in the tests of 
series II, the first speed listed is that at which the instability first 
occurred, but the other speeds are merely speeds above the first at which 
the frequency of the oscillation was measured. The last speed is that at 
which the tests of series II were discontinued; the nature of the insta- 
bility phenomenon did not appear to change in the speed range covered. 

The flutter-speed coefficient 2 v/Lwq,, the dimensionless dynamic 
pressure at flutter q*, and the frequency ratio f/fo; pertaining to the 
tests of series I are plotted in figure 3 against the density ratio p/pq- 

The frequencies measured in series II are shown plotted as a fimction of 
airspeed in figure i)-. In figure 5 are shown the speed and dynamic pres- 
sure (both in dimensionless form) at the onset of the yawing oscillations 
observed in the tests of series TV as functions of the density ratio p/Pq’ 

The flutter-speed coefficients 2v/ln^ and the frequency ratios f/fo; 

pertaining to the tests of series V are shown plotted as fuinctions of the 
density ratio P/Pq ia figure 6. 
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In figinre 7 the speed coefficients for which flutter or divergence 
occurred are shown plotted against Xq, for the tests of series VI. The 

points representing the various tests are not connected because the body 
mass and moment of inertia (and, hence, the still -air frequencies) were 
not constant in these tests. Figure 8 consists of a similar plot made 
for the tests of series VII. 

Observed flutter and divergence phenomena .- Several types of aero- 
elastic instability were encountered. In the tests of series I and V, as 
well as in some of those of series VI and VII, coupled flutter, similar 
to classical bending-torsion flutter was encountered, except that "bending" 
and "torsion" were lateral motion and yawing, respectively, in the case 
of these bodies. These two types of motion were distinguished by observing 
the strain-gage outputs . If the outputs of the front and rear gages had 
been in phase and of the same magnitude, the motion would have been from 
side , to side only, without yawing, but actually this type of motion was 
not observed in the tests. When the gage outputs were l80° out of phase 
and of the same magnitude, the motion was a pure yawing oscillation, and, 
when they were out of phase by any other angle, the motion was a coupled 
lateral-motion and yawing oscillation. 

In one test of series VII a combined flutter and divergence insta- 
bility was observed, not imlike the type of phenomenon which a wing may 
experience if its flutter occurs in a mode which involves very little 
bending; the tube tended to diverge to the stops after a few oscillations 
of increasing amplitude. As in all tests where divergence was observed, 
when the body began to diverge in yaw it moved over laterally as well as 
mder the action of the side forces brought into play by the yawed attitude. 

The flutter frequencies were ill-defined occasionally, particularly 
when the body at its flutter condition was also close to a divergence con- 
dition; that is, flutter then occiarred so suddenly that no definite fre- 
quency could be obtained from the strain-gage record. 

Observed uncoupled oscillations ..- An instability was observed in the 
tests of series II, III, and IV. This phenomenon consisted in continuous 
or intermittent, self -excited, small-amplitude yawing oscillations, usually 
with fairly well defined frequency. When the oscillations were intermit- 
tent they started up at random intervals rather than subsiding and increasing 
in a regular fashion, as do oscillations with beats. In two of the tests 
in which such oscillations occurred the frequency was insufficiently 
defined to be measured. This phenomenon differed from flutter not only 
in the fact that it involved small constant amplitudes but also in the 
fact that, unlike flutter, it involved no bending deformations of the 
struts and, hence, no lateral motions of the bodies. For lack of a better 
name this phenomenon is listed as "yawing oscillations" in table 3 and 
will be referred to as such hereinafter. If the body were mounted in such 
a way that the struts were horizontal, as would be the case if a body were 
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mounted on the side of a fuselage, this phenomenon would consist of 
angle -of -at tack or pitching oscillations. 

Some of the tests of series III which resulted in yawing oscillations 
were continued to speeds above that at which the oscillations started; the 
first speed listed is then the one at which the instability first occurred, 
and the second is that at which the tests were discontinued. Ko change 
was noted in the nature of the phenomenon within this range of speeds . 

The two frequencies listed for these tests correspond to these two. speeds; 
the values of p, q, and q* are those which correspond to the first 
speed. In the tests of series III in which no instability occurred, and 
also in the first test of series VT, the speed listed is that at which the 
tests were discontinued. 

Results of calculations .- Some of the results of. the instability cal- 
culations described in the appendix are also listed in table 3 and are 
shown in figinres 4, 6, 7> and 8. For series I and II the calculated 
flutter speed was infinite, that is, the calculations did not predict 
flutter for any finite speed. Nor was it possible to calculate a finite 
speed at which self -sustained yawing oscillations could exist, but the 
frequencies of the yawing oscillations of the body in a moving air stream 
could be calculated and are shown in table 3 and figure 4. 

Flutter speeds and frequencies were calculated for the tests of 
series V; the results are shown in table 3 and figure 6. Flutter speeds 
and frequencies as well as divergence speeds were also calculated for the 
tests of series VI and VII and are listed in table 3; the speed coeffi- 
cients are also shown in figures 7 and 8. 


DISCUSSION 

A Note Concerning the Speed and Dynamic -Pressure Parameters 


Two dimensionless parameters have been used in order to compare the 
results of the various tests. The first of these is substantially the one 
commonly used in flutter work, the flutter parameter or speed coeffi- 
cient v/bcjcj^, where b is the half -chord and is here replaced by one- 
half the (basic) body length, so that the parameter becomes 2 v/L6Dq,. It 

involves a measure of the dynamic pressure (in the term v, although the 
air density is not taken into account), as well as of the structmral and 
inertia properties involved in angular deformations (combined in the 
still-air yawing frequency (Dq,) . The other parameter used herein is one 

often used in static aeroelastic work and may be referred to as the 
divergence parameter or dimensionless dynamic pressure; the form of this 

parameter used in this paper is q* = ■ — - P . it represents the ratio of 
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the aerodynamic moment per unit angular displacement to the elastic 
restoring moment per unit displacement and Involves implicitly the dynamic 
pressure as well as the structural and aerodynamic properties pertinent 
to angular deformation, because the factor 2 is the value of the moment 
coefficient (referred to the vol\ame of the body) per unit angular dis- 
placement according to thin-body theory. This theory applies only to 
closed bodies, but, if it is extended to open bodies on the basis of the 
assimiption ths,t the flow inside the body has the same velocity as that 
outside and that the rear end of the body acts like a sharp trailing edge, 
the same value is obtained for the moment coefficient because the effects 
of these two assumptions cancel each other. 

The main advantage of the flutter parameter is that it includ.es some 
dynamic or inertia effects; on the other hand, the advantage of the 
divergence parameter is that by virtue of its explicit inclusion of aero- 
dynamic properties it serves as a more precise definition of certain 
instability phenomena. For instance, flutter and divergence can occur 
over a wide range of values of the flutter parameter, but, although flut- 
ter may occur at almost any value of the divergence parameter, divergence 
should occur at values of this parameter near unity, (if 2 is the correct 
value of the moment coefficient per unit angular displacement or if the 
correct value is used in the defininition of the divergence parameter 
instead of 2, divergence will occur when the parameter is 1.) Also, as 
shown in the appendix, the divergence parameter appears to play an impor- 
tant part in determining the frequency of the uncoupled yawing oscillations. 

Both parameters, therefore, have some advantages and, in view of the 
exploratory nature of the investigation, both have been used in attempts 
to analyze the observed phenomena. 


Flutter and Divergence 

Flutter and divergence similar to the classical instability phenomena 
on wings were both observed on the aerodynamically clean bodies, the 
streamlined body and the open tube with and without fairing on the struts, 
as may have been expected, because under certain conditions these bodies 
have linear aerodynamic forces which may be esqjected to give rise to 
phenomena similar to the classical instability phenomena of the wings. 

The streamlined body .- The streamlined body fluttered in the tests 
of series Ij in which the density was varied but all other parameters 
held constant, at speeds which corresponded to a wide range of the flutter 
parameter 2v/Lo^ but to values of q* which varied only between 0.79 

and 0.9^. The flutter frequency was substantially constant in these tests. 
(See fig. 5.) 
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The guantlty q.* varied relatively little in these tests because 
it contains the air density, which was the only variable in these tests. 

The small variation of q* in the tests can he considered to he due to 
the change in the mass ratio, that is, the ratio of the hody mass to the 
mass of the displaced air. Inasmuch as the hody would have diverged at 
q* = 1 if the actual moment on this hody in steady flow were that pre- 
dicted hy thin-hody theory, flutter occurred in these tests at dynamic 
pressures from 6 to 21 percent lower than the theoretical dynamic pres- 
sure at divergence. 

When attempts were made to calculate the flutter speed of the stream- 
lined hody, the aerodynamic forces predicted hy thin-hody theory were 
found to he incapable of predicting a finite flutter speed. For bodies 
with a fineness ratio of about 5^ "the moment due to angle of attack and 
normal force due to steady rotation predicted hy thin-hody theory are 
about 25 and 40 percent higher, respectively, than those predicted hy 
exact potential-flow theory. The values predicted hy potential-flow 
theory are, in t^lrn, somewhat higher than the actual values as the result 
of boundary- layer separation. Inasmuch as quantitative errors in the 
predicted forces would tend, hy themselves, to result only in an incorrect 
flutter speed, the fact that the predicted flutter speed does not even 
exist suggests that aerodynamic forces must he involved which are not 
predicted hy this theory. Such forces are the normal force due to angle 
of attack and the moment due to steady rotation. These forces are zero 
according to thin-hody and exact potential-flow theory, hut actually they 
do exist; often they vary linearly with angle of attack and rate of rota- 
tion, respectively, and thus represent the type of forces required to 
cause classical flutter. 

In the light of this discussion, prediction of the flutter speed of 
bodies of revolution thus requires a knowledge of the effects of the 
boundary layer and of the phenomena associated with its separation on 
the aerodynamic forces. Hence, the main shortcoming of thin-hody theory, 
insofar as the prediction of flutter is concerned, consists in the inability 
of potential-flow theory in general to predict some of the critical forces 
involved in these phenomena rather than in the degree to which thin-hody 
theory approximates exact potential-flow theory. In divergence, however, 
only the moment due to angle of attack is Involved; therefore, the correct 
magnitude of this force is important. Inasmuch as the moment predicted 
hy thin-hody theory is about 50 or ^4-0 percent higher than the actual 
value, the dynamic pressure at divergence would he that much hi^er than 
that estimated on the basis of thin-hody theory; therefore, in the tests 
of series I, the highest speed reached probably corresponds to about 
80 percent of the true divergence speed. 

The open tube on unfaired struts .- In the tests of series V the open 
tube fluttered at all densities. The values of 2v/Lcq^ decreased from 
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10.32 to 5 ’30 with increasing density, and the values of q.* appeared 
to increase gradually from O.516 to 0.579^ except for the second test of 
the series for which q* was somewhat hi^er than expected from the 
other three tests. (See table 3-) These trends are the same as those 
observed in the case of the streamlined body. In contrast to the stream- 
lined body, however, flutter speeds could be calculated for the open tube 
on the basis of thin-body theory (with the additional assumptions of 
unobstructed flow through the tube and of the flow continuing in the same 
direction after leaving the tail end of the tube rather than reassuming 
the free-stream direction) , These speeds are in fairly good agreement 
with the measured ones. 

The observed frequencies did not vary with the density, a fact which 
was also noted in the case of the streamlined body. The calculated fre- 
quencies did not vary with density either but were about 30 percent higher 
than the observed ones. 

On the basis of these comparisons it appears that the aerodynamic 
forces are taken into account in a qualitatively correct manner, but that 
quantitatively they must be improved considerably before they can be used 
to predict flutter speeds and frequencies correctly for ducted bodies. 

In the tests of series VI, the tube fluttered or diverged in all 
cases, except in the first test, in which the speed of the test was not 
carried to a high enough value. Flutter occurred when the center of 
gravity was at or behind the elastic axis, and divergence, when it was 
at or in front of the elastic axis (see fig. ’J) } when it was at the 
elastic axis, the tube fluttered at the higher mass and diverged at the 
lower mass. This trend agrees with the trends noted in the tests of 
series I and V, because in these tests there appeared to be a tendency 
to approach divergence as the air density increased and, hence, the mass 
parameter decreased. The values of 2v/lu(^ at instability varied from 

5.^7 to 7-95 and those of q* between 0.459 and 0.754 (see table 3). 

The values of q* at divergence tended to be higher than those at flut- 
ter, as may be expected in view of the nature of q* as, primarily, a 
divergence criterion. No such distinction can be made in the case of 
the values of 2v^In^ at Instability, both the highest and lowest values 

of which correspond to divergence. 

The calculated speeds at instability are in fairly good agreement 
with the measured speeds for the rearward locations of the center of 
gravity, but the flutter speeds predicted for forward center -of -gravity 
locations are much too low. The two measured flutter frequencies are 
substantially below the values calculated for those two cases. The cal- 
culated divergence speeds are consistently higher than the measiired ones 
by about 20 percent on the average; therefore the moment coefficient per 
unit angular displacement (the only aerodynamic parameter which enters 
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the calculation of the divergence speed) must he about 40 percent higher 
than that estimated by means of the modified thin-body theory used herein. 
Actually the tube tinder consideration is not sufficiently slender to 
justify the use of thin-body theory, and its abrupt changes in cross- 
section at the nose and tail violate certain assumptions Inherent in 
thin-body theory; also, the validity of the two additional assumptions 
is doubtful. Hence, the inability of the modified thin-body theory to 
predict q^uantitatively useful results is not surprising. 

A more accurate potential-flow solution for the flow through and 
about an open tube would be difficult to obtain, and its validity would 
still be open to question because the flow inside the tube would be 
significantly affected by the boundary layer in the inside walls. The 
available solution for a ring airfoil is inapplicable to this case 
because the tube is far too slender for this theory. The most promising 
solution, therefore, appears to be the use of a semiempirical method for 
estimating the aerodynamic forces required in flutter analysis. Such a 
method might consist in retaining thin-body theory but modifying the two 
additional assumptions, that is, by estimating the magnitude of the 
forces which are, in effect, concentrated at the rear of the afterbody 
and the extent to which the flow decelerates inside the tubes on the 
basis of measurements of the moment and normal force due to angular dis- 
placement. In divergence calculations only the moment per unit angular 
displacement is required, of course. 

The open tube on faired struts .- The aeroelastic instability phenom- 
ena of the tube were substantially unchanged by the addition of the 
fairing; apparently, the increase in the aerodynamic forces was canceled 
in effect by the increase in the stiffness of the configuration. Flutter 
still tended to occur at the further rearward position of the center of 
gravity, and divergence, at the forward positions. (See table 5 and 
fig. 8.) The agreement between the calculated and measured speeds was 
poor, the calculated values being much too low, and the calculated fre- 
quency corresponding to the one measiired frequency was also much too low. 
The first two tests of this series serve to illustrate the difficulty of 
estimating the speed at which aeroelastlc-instabillty phenomena occur; 
under identical test conditions, the model diverged in one case at 
522 feet per second and exhibited some symptoms of flutter at 48l feet 
per second in the other case, flutter apparently having been suppressed 
the first time. 

The fact that the calculated values of the divergence speed are 
lower than the observed values indicates that, inasmuch as the aerodynamic 
moment on the tube due to angular displacement was probably imderesti- 
mated, as in the tests of series VI, the forces, on the fairing were 
greatly overestimated, as might be expected to be the case because unmodi- 
fied two-dimensional theory was used to estimate these forces. The exact 
calculation of the mutual interference effects of the tube and the wing 
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tip represented by the fairing represents a very difficult problem for 
which no solution has been foxmd to date. For the time being, therefore, 
approximations similar to those made herein must be- used, although they 
could probably be improved by resorting to empiricisms based on some 
measured results. 


Yawing Oscillations 

The streamlined body; speed at which oscillations start .- When in 
the tests of series II the streamlined body was tested on stiffen struts 
and with an elastic -axis location 5 percent of the body length further 
forward than in the tests of series 1, no flutter was observed; instead, 
the body experienced the self -excited yawing oscillations described pre- 
viously. These oscillations started at a speed which is relatively low 
compared to those at which flutter occurred in the tests of series I; 

2v 

it corresponds to - — = 3.76 and q.* = 0.211, and the speed at which 

I<u 

a 

the tests were discontinued corresponds to = 6.70 and q.* = O. 63 O, 

whereas at the same density flutter would have occurred in the tests of 
series I at values of 2vy^Lcu^ and q* of about 8.0 and 0.9^ respec- 
tively. Therefore, if the tests of series II had been continued to a 
speed some I 5 or 20 percent higher than that at which they were discon- 
tinued, flutter would probably have occurred if the values of 2v/Lo)^ 
and q* at flutter are assumed not to differ much between the two test 
series. The fact that oscillations occurred in the tests of series II 
but not in those of series I is probably the result of the difference in 
elastic -axis locations, as will be shown later. 

In order to calculate the speed at which yawing oscillations might 
start, an attempt was made to solve the equations of motion for the case 
of zero lateral displacement. For this case there are two differential 
equations with one imknown function. W'ith the aerodynamic forces given 
by thin-body theory, however, the equations can have a solution only when 
the airspeed is zero. Therefore, the nature of the air forces must again 
be different from that assumed to yield equations which are compatible 
at nonzero airspeed. 

From a physical point of view, self-excited oscillations can occur 
only if the net damping is zero. In the equation for the equilibrium 
of the moments as set up in the appendix, there is no damping term, 
because the aerodynamic damping moment according to thin-body theory is 
zero and the structural damping was assumed to be zero. If, however, 
the aerodynamic moment were actually small and negative, at a certain 
speed it would be just large enough to counteract the structural damping, 
and self -excited oscillations would start at this speed. At higher 
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speeds^ the net damping would he negative, hut nonlinearities in the 
aerodynamic forces might prevent the oscillations from diverging. Inas- 
much as these aerodynamic forces would then he due entirely to deviation 
from potential flow, they would prohahly vary with Reynolds number so 
that the q^uestlon whether and when a hody might experience self -excited 
oscillations prohahly depends on the Reynolds nimiher involved, as well 
as on the structural damping. 

The streamlined hody; frequency of oscillations .- Although the equa- 
tions set up for the analysis of the yawing oscillations do not furnish 
a solution for the speed at the onset of the oscillations they do give 
an indication of the frequency of the resulting oscillations, that is, 
of the "natural" frequency in moving air. This frequency is expressible 
as a product of the still-air yawing frequency, which involves the dynamic 
characteristics of the system, and a correction factor which, except for 
a generally negligible dependence on the mass ratio, is a function only 
of the static aeroelastic characteristics represented in the param- 
eter q*, namely \Jl - q*. That the trend of the calculated frequencies 
agrees well with the trend of the frequencies measured for the stream- 
lined hody may he seen from figure A. However, the rate of decrease 
with speed is less than predicted. This fact suggests that, inasmuch 
as the expression for the frequency is obtained directly from the equa- 
tion representing the equilibrium of the moments on the hody, the esti- 
mated aerodynamic moments on the streamlined hody are too large, which 
is true, as previously noted. 

Miscellaneous hlirff bodies; effect of elastic -axis location on the 
speed at which oscillations start .- Yawing oscillations occurred in the 
majority of the tests of series III. The speeds at which they started 
correspond to values of 2v/Lo:^ and q* much lower than those at which 

the streamlined hody fluttered in the tests of series I and of about the 
same magnitude as those at which the oscillations of the streamlined hody 
started in series II. (See table 3.) Some of the tests of series III, 
for instance, the last-listed one, resulted in no instability at what 
appear to he fairly large values of 2v/ln^ and q*; however, all values 

of 2 v^Lo)q, listed in table 3 are based on the length of the basic tube. 

If they had been based on the actual lengths of the various bodies, 

2vy'Lo^ would have been about for the last-listed test, which value is 
about the same as that at which the oscillations of the streamlined body 
began (3-76) and much lower than that at which that body could have flut- 
tered at this density (about 8). Also, although the value of q* (O.76O) 
at which this test was discontinued is high compared with that at which 
the oscillations of the streamlined body started (0.211), it is lower 
than that at which that body would probably have fluttered at this density 
(about 0.9) • The values of q* in table 3 for series III are based on 
the actual volume, as are those for the other series. 
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The speeds at which oscillations started, in the form of the parame- 
ter 2v/ln^, are shown plotted in figure 9(a) as functions of the dimen- 
sionless elastic-a^is location a. The importance of the elastic-axis 
location in determining the speed at which the oscillations start may be 
deduced from the equation set up in the appendix for the equilibrium of 
moments in the 3 rawing-osc illations phenomenon, which contains the dimen- 
sionless elastic-axis location a. 

In figure 9(a) the upward -pointing arrows on several sketches refer 
to cases in which no instability was observed; therefore, any instability 
would have had to occur at values of 2vjlxa^ greater than those shown, 

which represent those at which these tests were discontinued. Figure 9(a) 
also contains one point representing the conditions at which the stream- 
lined body started to oscillate in the tests of series II. The sketches 
in figure 9(a) all imply an airstream approaching from the left. 

The only body for which sufficient information was available to 
deduce the effect of elastic-axis location on the speed at which the 
oscillations started is the one with hemispherical nose and tail. In 
figure 9(a) this body is represented by four points which appear to lie 
on a smooth curve which has a minimum at a = 0. For the body with hemi- 
spherical nose and streamlined tail, two points were available. These 
points are connected by a curve based on the pattern exhibited by the 
body with hemispherical nose and tail which represents the probable 
variation of the speed coefficient with elastic-axis location, although 
the minimum of the curve may not be at a = 0 as shown. For each of 
three other bodies (including the body used in the tests of series I 
and II), one point was available which represented the onset of oscilla- 
tions and one point which indicated only that the oscillations, if any, 
would have to start at higher speeds except that, for the steamlined 
body, the second point represented the flutter condition. Curves repre- 
senting the estimated variation of speed at the onset of oscillations 
with elastic -axis location are shown for two of these bodies as well. 

These approximate curves indicate that the speed at which oscillations 
start is lowest when the elastic axis is near the midpoint of the body. 

^ using the estimated variations as a guide, that is, by estimating 
what the speed would have been if the elastic-eixis -location parameter a 
had been zero, the effect of the body shape on the speed at which oscil- 
lation's tend to start can be divorced from that of the elastic -axis loca- 
tion. 


Miscellaneous bluff bodies; effect of body shape on the speed at 
which oscillations start .- Inasmuch as a part of the large differences 
in the speed coefficients shown in figure 9(a) Tor the various bodies is 
due to the fact that the coefficients were based on the length of the 
basic tube, the speed coefficient will be considered to be based on the 
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actual length of the body in the following discussion. The various 
bodies then fall into several main classes. Oscillations appear to start 
at the lowest speeds when the tall of the body is hemispherical, regard- 
less of whether the nose is hemispherical, square f cover plate), or a 
small cone (for all of which is about 1.5 or a streamlined 

body (for which 2v is about 1.5)» The class with the next highest 

range of speeds at which oscillations start is that with tails consisting 
of the streamlined shape or the large cone. If for this group the nose 
is hemispherical or square, 2Yjlxa^ is about 2.5; if the entire body is 

streamlined, 2 v/LxIq, is about 3; and if the nose of the body is a pointed 
streamlined shape or a large cone, 2 v/Lo^ is greater than 4. The class 

for which oscillations start at the highest speeds is, surprisingly, that 
with a square tail (cover plate); if for this group the nose is hemispheri- 
cal or square, 2v^Inj^ is about 5 to 4, and, if the nose is a pointed 

streamlined shape, 2 v^Lo)q, is greater than 4. All these nmerical values 

must be used with caution, of course. 


In general, then, the speed at which oscillations start appears to 
depend to a large extent on the shape of the tail of the body, a hemi- 
spherical tail being the least favorable in that it oscillates at the 
lowest speeds, a streamlined tail or large cone used as a tail being much 
better, and a square tail being most favorable in this respect. The 
shape of the nose has almost no effect when the tail is hemispherical but 
has some effect in the other cases, a square or hemispherical nose being 
worst, a conventionally streamlined (rounded) nose being better, and a 
pointed streamlined nose, best. The aerodynamically cleanest configura- 
tion, the streamlined body, occupies a relatively favorable place; the 
speed at which it may start oscillating can be increased further by 
replacing the rounded streamlined nose by a pointed one. 

The only way in which the numerical values given for 2 v/Lcd^ in the 

preceding discussion can be related to those corresponding to the classi- 
cal Instability phenomena is by noting that the streamlined body fluttered 
at values of 2Yjlm^ of about 8. As a rule of thumb, then, based on these 

very limited data, an aerodynamically clean body may be expected to start 
oscillating at speeds as low as about one-third its flutter or divergence 
speed. To attempt a similar correlation for bluff bodies would be futile, 
because these bodies do not experience the linear aerodynamic forces on 
which the classical instability phenomena are premised. 


The reasons for the relative behavior of the various bodies are as 
yet somewhat obscure. The effects of the nose on the speed at which 
oscillations start are probably associated with separation at the nose, 
because the less disturbance caused by the nose, the more favorable the 
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cdnf iguration. Similarly, the behavior of the tail can probably be 
explained by the effects of separation at or ahead of the tail. The 
unfavorable effects of the hemispherical tail, for instance, are probably 
the results of the rather large projected area of the tail (projected on 
a plane throu^ the axis of the body), which is exposed to the separated 
flow. The relatively favorable sq.uare tail, on the other hand, has no 
projected area exposed to this flow. The streamlined tail probably 
causes relatively little separation and thus is relatively favorable if 
the nose is favorably shaped; this effect corresponds to the observed 
effects . 

In view of the complicated nature of these phenomena there appears 
to be little hope of arriving at theoretical methods of predicting the 
speed at which the oscillations start, at least in the case of bluff 
bodies which are unlikely to be exposed to high-speed airstreams anyway; 
in the case of more or less streamlined bodies, empirical approaches may 
prove fruitful. In either case the oscillations do not appear to be vio- 
lent, and, if they are undesirable, they can always be eliminated by 
stiffening the supports and in many cases merely by changing the body 
shape or the elastic -axis location, and possibly also by using vanes 
mounted on the body. 

Miscellaneous bluff bodies; frequency of oscillations .- The ratio 
of the freq.uency of the oscillations to those in still air is shown 
plotted in figure 9(ti) as a function of the dimensionless dynamic pres- 
sure q.* corresponding to the speed at which the oscillations started. 

The arrows in figvire 9(^3) refer to cases in which no frequency was meas- 
ured and serve merely to call attention to the fact that oscillations 
did start on the given model at the indicated value of q*. As in fig- 
ure 9{a.) , one point in figure 9(ti) represents the conditions at which 
the streamlined body started to oscillate in the tests of series II. 

Also, as in figirre ^{a.) , the sketches in figure 9(l3) imply an airstream 
approaching from the left. 

The frequencies shown in figure 9(^) agree fairly well with the theo- 
retical curve (which neglects the effect of the mass ratio) . In view of 
the fact that many of the assumptions made in the analysis are violated 
by the bluff bodies, this agreement is better than may have been expected. 
All points in figure 9(b) pertain to tests at densities close to sea 
level. The results of the tests of series IV at various densities follow 
a similar pattern, although the range of q* covered by these tests is 
small. (See fig. 5.) The frequency of the oscillations can thus be 

estimated on the basis of the relation — = \jl - q* with fairly good 

% 

accuracy. If the body is aerodynamical ly clea,n and the speed relatively 
high, the accuracy of this formula can be improved by replacing the fac- 
tor of 2 in the definition of q* by a better value of the moment coef- 
ficient per unit angular displacement. At low speeds, the frequency may 
be expected to be substantially the same as the still -air frequency. 
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CONCLUSIONS 


1. Streamlined "bodies and open tubes mounted on thin flexible struts 
which do not contribute any aerodynamic forces have been found to diverge 
and flutter I flutter tended to occur for relatively far rearward locations 
of the center of gravity and for relatively high mass ratios (body mass to 
mass of displaced fluid) . 

2. Flutter could not be predicted for the streamlined body by using 
aerodynamic forces based on potential theory. For the open tube, the 
calculated flutter and divergence speeds did not agree well with the 
measixred values; the discrepancies are believed to be due to the intrinsic 
shortcomings of potential -flow theory. The analysis of imsteady aero- 
elastic effects of bodies of revolution therefore appears to req.uire a 
knowledge of the boundary- layer and separation effects on the unsteady 
forces on these bodies. 

5. Closed strut-mounted bodies of revolution appear to be subject 
to a nonclassical instability which consists in self -excited nonviolent 
oscillations which appear to start, in the case of aerodynamically clean 
bodies at least, at speeds about one-third that at flutter or divergence 
for the given body, 

4. The speed at which the oscillations start for a given body depends 
on the elastic-axis location, being lowest when the elastic axis is 
located near the midpoint of the body. This speed is also determined to 

a large extent by the shape of the body, particularly of the tail end. 

For bodies with hemispherical tails, the oscillations start at low speeds 
but, for bodies with streamlined and, particularly, with sq.uarely cut off 
tails, they start at relatively high speeds; the optimum nose shape appears 
to be a pointed streamlined shape and the next best, a conventional stream- 
lined nose shape. 

5. The mechanism which causes these oscillations is as yet unknown, 
although negative aerodynamic damping moments appear to be req,uired. The 
speed at which the oscillations start cannot, therefore, be predicted at 
present; its calculation apparently requires a knowledge of boundary- 
layer and separation effects on the unsteady aerodynamic forces. The 
frequency, however, can be estimated from a simple formula Involving the 
frequency of the oscillations in still air and the ratio of the given 
dynamic pressure to that at divergence. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May 12, 1955. 
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APPENDIX 

CALCULATION OF AEROELASTIC INSTABILITIES 


Equations of Motion 


The equations of motion of a rigid body mounted on flexible supports 
and performing combined yawing and lateral oscillations are 


-iQa - m ^2 ” ®i^^ - (l + ig)K^^ + = 


0 


( 1 ) 


-m^S 2 - - mil - (l + ig)Kj^h + P = 0 


( 2 ) 


where the dots designate differentiation with respect to time, and where 
P is the aerodynamic force and the aerodynamic moment (positive in 

the direction of positive h and a, respectively) . The manner in which 
P and were calculated for the bodies considered in this paper is 

described in the following sections. 


Aerodynamic Forces on Closed Bodies 

The aerodynamic forces on closed bodies of revolution performing 
unsteady motions in supersonic flow have been calculated by linearized 
potential-flow theory (see refs. 5 and 6). The aerodynamic forces on 
bodies of revolution in steady incompressible flow can easily be calcu- 
lated by potential -flow theory (by using sources and sinks, for instance, 
as in ref. 7)j according to this theory, the normal force is zero. The 
exact calculation of the aerodynamic force and moment for unsteady motion 
by potential -flow theory, however, is quite difficult both at subsonic 
and supersonic speeds, and in view of the fact that they are known to be 
influenced to a large extent by the effects of viscosity, a large expendi- 
ture of effort in calculating them is hardly warranted. In the absence 
of any means of taking the effects of viscosity into account for unsteady 
motions, a simple approximation to potential-flow theory, namely slender- 
body theory, has therefore been used for the purpose of the calculations 
described herein. (See ref. 8, for instance, for an outline of a slender- 
body theory in quasi-steady flow.) 

The assumption made in slender -body theory is that the momentum of 
the flow in a plane perpendicular to the free stream is the same as it 
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would be if this flow were two-dimensional. This assumption implies 
that the derivative of the radius with respect to distance along the 
length of the body is small (which implies, in turn, that the body is 
very slender, that is, the body has a high fineness ratio), and, also, 
that the angle of attack and any motions are small. (For a fuller dis- 
cussion of these assumptions from the mathematical point of view, see 
ref. 9.) 

The momentimi of the flow about a circular cylinder for a unit length 
along the cylinder is equal to the product of the rate of motion of the 
cylinder and the apparent mass, which is equal to the mass of air dis- 
placed by a unit length of the cylinder, or pnR^. At any section of a 
body of revolution the rate of motion relative to the component of the 
free-stream velocity normal to the axis of the body is va + fi + (s - Sp)o: 

this value is within the approximations implied in slender -body theory. 
Therefore, the momentiam per unit length along the body of the flow per- 
pendicular to the body at this section is pjtR2(s) jva + h + ^s - s-j^j aj . 

The force exerted by the body on the fluid per unit length along the 
body is equal to the time rate of change of the momentum per unit length, 
the rate of change being that along the path of a particle, that is. 

But within the approximation Implied in slender-body theory, 

= 1. + V-^ 

Dt dt Sx 


Therefore, the force per unit length along the body is 


1 ., 2 / ■ r 


I “ V— j ■jprtR'^(s) l^va + h + ^s - 


Hence, if the body is performing oscillations defined by 


iojt 


a = a^e 


h = h^e 


iojt 


( 3 ) 

(h) 

(5) 
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then 


r 


I = -p ( io) + V itR^ ( s ) 

os 


va + iaho + ico^s - s^^oIq 


i(X!t 


-P. 


2io3v«R^(s) - 0 ^ fs - s^]rtR^(s) + — itR^(s) + 


ivoi^s - itR^(s) 


«o + 


-a:i%^(s):t + ivcu— jtR^(s) 

Bs 


o 


^ iojt 
> e (6) 


The aerodynamic force P (positive in the same direction as h) 
and moment (positive in the same direction as a) are then 


P = / Z ds 

Jo 


= -pVk ^ 


iojv + CD 


'K - 5) 


- ca^h 
O r 


ioyt 


(7) 


and 


% = j (s - Sj^)l ds 


PV^ <: 


V^ + CD^(s" 


2sj^s + S]^ )|q;^ + 


CD 


■(s - S^) + i 


IVCD 


ho!>e 


icDt 


(8) 


vhere 


s = -i- / sjtR^(s) ds 

^h Jq 
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and 


= “ r s^jtR^(s) ds 


''Td-'O 


or, in dimensionless form. 


py2L2gioJt 


2iklQ + 4k2(alQ - Ii) 


2^ 2ho 


( 9 ) 


and 


M 


a 


pv 


2L3ei“fc 


Iq + - 2al^ + I 2 ) 


% + 


ikio + 2k^(li - oIq) 


2h. 


( 10 ) 


where Iq = / S(|) d|; I = / |S(|) d|; and Ip = / I S(|) d| . 

^0 O 0 Jo 


For the airfoil-shaped body of revolution, these three values are 
Iq = 0 . 01626 , Ip = 0 . 00674 , and I 2 = 0.00535. 

For steady flow (o) = 0) equations ( 7 ) and (8) give the known results 
of slender -body theory 


P = 0 


( 11 ) 


= pv^v^a 


( 12 ) 


Aerodynamic Forces on Open Bodies 

In attempting to calculate the forces on the open tube in a similar 
manner, several problems arise. If the recovery factor is assumed to be 
100 percent, that is, if the velocity of the flow in the tube is assumed 
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to be eqixal in magnitude to the free-stream velocity, the combined appar- 
ent mass of both the internal and the external flow at a given section 
of the body is just twice the apparent mass of the external flow used 
heretofore. However, one of the assumptions of slender -body theory, 
namely, the one concerning continuity of the radius along the length of 
the body, is violated at the nose and tail section of the tube. Conse- 
quently, as a result of the abrupt changes in cross-sectional area, not 
only the concentrated forces predicted by slender -body theory at the nose 
and tail section but also the distributed forces predicted on the remain- 
der of the tube are open to question. 

For lack of a better theory a modified slender-body theory 1ms been 
used to calculate aerodynamic forces on the open tube. The modification 
consisted in disregarding the concentrated forces on the tail section on 
the premise that both the external and the internal flows leave the 
trailing edge of the tube tangentially and are not realined with the free 
stream. This assumption is equivalent to the K\itta condition of subsonic 
airfoil theory and is used also in the application of slender-body approxi- 
mation to airfoil theory. (See ref. 10.) The assumption is thus, essen- 
tially, that the exit section of the tube acts like the trailing edge of a 
wing of very low aspect ratio. 

With these approximations equations (6), (9)^ (lO), (ll), and (l2) 
become, in the case of an open cylindrical tube. 



where 6(s) is the delta function which represents the concentrated 
loads , and 
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^ — = - 2 rt J 21k(l - 0 ^) - a - k-^(- - a + 


a + 
o 




-2k^(i - a) - ika 



L 


w 


(15) 


^o>=0 = -2pJtR^v^a 


(16) 


M = 2prtR^s-, v^a 


(IT) 


Calculation of Flutter Speeds and Frequencies 

For the streamlined body, the force and moment coefficients given 
by equations (9) and (lO) vere introduced into equations (l) and (2) 
with a and h given by equations (4) and (5), and an attempt was made 
to solve these equations by the conventional methods of two-degree-of- 
freedom flutter analysis. (See ref. k, for Instance.) However, no 
solution was foimd to exist; therefore, if the aerodynamic forces were 
correct, the body would not experience flutter at any finite speed. For 
the open tube the force and moment coefficients given by equations (l4) 
and (15) were substituted into equations (l) and (2). In this case 
flutter speeds and frequencies did exist and the computed values are 
given in table 3 and are shown in figures 6 to 8. 

For the tube on the faired struts the forces and moments given by 
equations (l4) and (15) were used for the tube proper. For the fairing 
the force and moment were assimied to be given by two-dimensional theory 
at any section and were obtained from reference 4. The aerodynamic inter- 
action between tube and fairing was thus neglected, as was the effect of 
the finite span on the forces experienced by the fairing. A Raylelgh-Ritz 
type of analysis was used with two modes, a linear torsion and a parabolic 
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"bending mode. These modes were selected on the basis of the consider- 
ation that the struts deflect somewhat as shown in the two following 
schematic front views: 



Lateral displacement 



The flutter speeds and frequencies calculated in this manner are given 
in table 3 and are also shown in figure 8. 


Calculation of the Frequency of the Uncoupled Yawing Oscillations 


In order to determine what characteristics, if any, of the yawing- 
oscillation type of instability could be predicted, cognizance was taken 
of the fact that these oscillations did not involve any bending deflec- 
tions; hence, h was set equal to zero in equations (l), (2), (9), 
and (lO) and, for the sake of convenience, the structural damping coef- 
ficient was assumed to be zero as well. The following equations resulted 




2t , 2^.3 

> + pv L 

” “I 

Iq + 4k^(o^l0 - 2al^ + I 2 ] 


w/ 


(18) 
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and 


r- 




2t2 

pv L 


2iba;o + 


’(“lo - 



A solution of equation (l8) can be written as 


(19) 



( 20 ) 


This solution furnishes no information concerning a speed at which the 
oscillations may start but, instead, gives the frequency at which the 
body will tend to oscillate if it is yawed away from its equilibrium 
position and then released. This frequency is also the dominant fre- 
quency of the response of the body to random excitation. The term 


pl5/ 2_ 

la I 

I \ 
a 


0 


- 2al^ 



in the denominator of the expression is the ratio 


of the moment of inertia of air at free-stream density contained within 
the body to the moment of inertia of the body alone, both taken about the 
elastic axis. This term is inversely proportional to the mass ratio; it 
depends on the shape of the body to some extent but is substantially 
independent of the elastic-axis location. Except at extremely low mass 
ratios, this term is negligibly small; for the airfoil-shaped body for 

instance, it is 0.0020 -2- and 0.0019 ~ for the elastic -axis location 

Po Po 

used in the tests of series I and II, respectively. 


The numerator of the expression under the radical in equation (20) 
is equal to 1 - q*, where q* is defined by 


q* = 




( 21 ) 
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Therefore, if the small term in the denominator is neglected, eq.ua- 
tion (20) can he written as 


f 

— = \jl-q* (22) 

a 


In order to determine the speed at which these oscillations should 
occur equation (l8) must be solved simultaneously with equation (I9) • 
However, the only real solution of equation (19) is V = 0 and, in 


addition. 


2pL3 



I^ j. Therefore, if the aerodynamic forces 


given by equations (9) and (lO) were correct, oscillations could occur 
only at zero airspeed. These oscillations would then be the ordinary 
still-air yawing oscillations, the condition on Xq, being that necessary 


to uncouple the yawing from the sidewise-bending mode. 


Calculation of Divergence Speed 


Inasmuch as divergence is a static instability phenomenon, the speed 
required to diverge the body can be found by setting h, to, and q 
equal to zero in equations (l) and (lO) or, more simply, setting o) = 0 
in equation (I8) . Thus, 


or 


^divergence 




(23) 


divergence 


K 


= -2L 


The parameter q* defined in equation (21) is thus equal to the ratio 
of q to ^divergences dynamic pressure at divergence being that 


calculated by using slender -body theory. 


The divergence speeds for the tests of series VI and VII are given 
in table 3 and are shown in figures 7 and 8. For series I and II the 
value of ^divergence obtained from equation (24), and for any 
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other test it can he obtained fran the values of q* and q given in 
table 3. Therefore, 


‘^divergence = 1*97 Kq, for series I and II 
= 29-^ for series I 
= for series II 

= in general 
q* 
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TABLE 1.- ORDINATES USED TO GENERATE 
THE AIRFOIL-SHAPED BODY 


s/L 

r/l 

0.005 

0.0134 

.0075 

.0186 

.0125 

.0237 

.025 

.0315 

.05 

.0435 

.075 

.0529 

.10 

.0608 

.15 

.0732 

.20 

.0829 

.25 

.0901 

.30 

.0954 

.35 

.0985 


.0999 

.45 

.0993 

.50 

.0960 

.55 

.0906 

.60 

.0829 

.65 

.0736 

.70 

.0631 

•75 

.0516 

.80 

.0397 

.85 

.0277 

.90 

.0162 

.95 

.0061 

1.00 

0 
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TABLE 2.- DIAMETERS AND EFFECTIVE SPRING 
CONSTANTS OF STRUTS 


Strut 

Strut diameter^ 
in. 

Kh 

% 

A 

0.014-0 

10 

20 

B 

.100 

50 

250 

C 

.166 

110 

I4-5O 

D 

.251 

200 

820 
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TABLE 5.- RESULTS OF AEROELASTIC-INSTABHITY TESTS 


Bocly 

Nose Tail 

Strut 

m 


a 



fa 

Type of 

instability P X ICP 

V 

exp 



f 

e3(5> 

f 

fttlAl 


(^)anal 

* 

exp 

■ 


— 

r* 1 




Test series I 










A 0.1501 0.0690 -0.28 0.10 1.20 


B 0.1500 'O.O67O -0.18 0 4.00 



0.1578 

.1578 

.1555 

0.1014 

.1034 

.1091 

-0.44 

.24 

-.20 

0.44 

-.24 

.20' 

3.08 

4.03 

5.85 

WESm 

.1264 

0 

.06 


.1274 

<^.20 

.26 

5.61 

.1635 

.1635 

.1626 

.1910 

0 

.44 

.24 

-.20 

5.96 

4.11 

.1392 

.1592 

.0904 

.1045 

0 

-.20 

-.04- 

.16 

4.54 

5.85 

.1365 

.0926 

-.20 

.20 

5.91 

.1394 

.0950 

-.20 

.16 

5.86 

.1548 

.0926 

-.20 

.20 

5.87 

.1677 

.1630 

.28 

0 

5.92 

.1583 

.0940 

-.20 

.20 

5.90 

.1436 

.1180 

-.20 

-50 

5.65 

.1677 

.1738 

0 

-.26 

4.00 

.1635 

.1766 

0 

-.24 

3-92 

.1776 

.2100 

0 

-.18 

5.16 


Test series II 


Yaving oscillation 


Test series III 


7.9 Yaving oscillation 

0.2 Yawing oscillation 

4.4 Yaving oscillation 

9.2 Yaving oscillation 

2.9 None 

8.0 Yaving oscillation^ 

.5 Yaving OBClUatlon^ 

.9 Yaving oscillation^ 

.4 None 


None 

Yaving oscillation 


None 2 .o 4 524 

Yaving oscillation'’ 2.31 211-301 

Yaving oscillation 2.27 225-256 

Hone 2.21 368 

None 2 .18 424 

.8 Yaving oscillation 2.28 156-236 

.2 None 2.15 479 


31.4 1.70 

• 15.4 

31.8 1.75 

12.5 

31.2 1.71 

10.6 

32.2 1.69 

9.5 

33.2 1.60 

8.4 

37-0 (a) 

7.8 



104 ■ 


155 


171 


222 


282 


510 


7.5 

10.3-10.1 

14.1-13.9 

8.58-8.34 


7.8 

14.7-14.5 


3.84 

1.95-2.91 

1.97-2.58 

3.25-'t.l6 

4.75 

3.90-5.16 

7.39 

3.65-4.82 

‘^.75 

5.45 

1.93-2.22 

4.48 

3.32-4.73 

1.92-2.20 

5.55 

6.75 

2.55-4.18 

8.48 


^ot obtained. 
'’Intermittent . 
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tabu: 5.- KESULTS OF AEROELASTIC -INSTABILITY TESTS - Concluded 


Body 

Nose Tail 

Strut 

m 

la 

i 

a : 

1 

1 

“*a 

■ 

fa 

Type of 
instability 

p X lo5 

^exp 

^anal 


f 

■^exp 

^anal 

(^)exp 

'^/anal 

* 

^ exp 







Test series 

IV 











0.25 

423 



20.6 

10.5 


5.07 


.k2 

314 

— 

20.7 

10.7 


3.78 


.61 

253 

— 

16.6 

10.5 


2.80 


.80 

194 

— 

15.1 

10.7 


2.34 


.99 

171 

— 

14.5 

10.4 


2.05 


1.19 

160 

— 

15.2 

10.7 


1.92 


1.38 

1-57 

160 

162 

— 

17.7 

20.6 

10.7 

10.7 


1.92 

1.94 


1.77 

152 

— 

20.5 

10.6 


1.82 


1.96 

158 

— 

24.4 

10.6 


1.90 


2.15 

145 

— 

22.6 

10.6 


1.74 


2.31 

129 

— 

19.2 

10.4 


1-35 


0.1378 


0.1501 


4.35 


10.6 


Yawing oscillations 


0.053 

.054 

.043 

.059 

.037 

.039 

.046 

.053 

.053 

.064 

.058 

.050 


Test series V 


0.1410 


0.0864 


0.14 


1-37 


2.40 


Flutter 



r 

0.55 

195 

l 4 o 

10.5 

1.51 

1.99 

10.32 

7.42 

0.516 



1.12 

147 

103 

12.1 

1.54 

2.00 

7.80 

5-50 




1.55 

117 

88 

10.6 

1.58 

1.96 

6.19 

4.68 

.521 



2.35 

100 

71 

11.8 

1.56 

1.99 

5-30 

3.76 

.579 


Test series VI 


0.1237 

0.0681 

0 

0 


10.1 

None 

2.31 

323 



121 



2.04 


0.471 

.1410 

.0875 

0 

.08 

5»io 

8.4 

Flutter 

2.14 

377 

“296 

152 

(a) 

6.8 

5 .71 

4.48 

.596 

.1410 

.0708 

0 

0 

3-11 

9.2 

Divergence 

2 .l 4 

396 

=312 

168 

^■■3 

5.47 

6.75 

.659 

.1410 

• 0957 

0 

.Ik 

3.07 

8.2 

Flutter 

2.14 

369 

“304 

146 

5.50 

6.7 

5-73 

4.72 

.573 

.1580 

.1217 

0 

.26 

3.01 

7.4 

Flutter 

2.14 

353 

°323 

133 

(a) 

6.1 

6,06 

5.56 

.522 

■ 1750 

.1475 

0 

.36 

2.81 

7.0 

Flutter 

2.14 

331 

=351 

117 

4.60 

5.6 

6.00 

6.38 

.459 

.1580 

.1217 

0 

0 

5.00 

7-3 

Flutter 

2.14 

389 

“276 

162 

(a) 

6.2 

6.78 

4.82 

.636 

.1750 

.1475 

0 

-.14 

2.83 

6.8 

Divergence 

2 .l 4 

424 

=177 

192 


^ 6.6 

7-95 

9.13 

.754 

.1580 

.1217 

0 

-.26 

2.96 

7.5 

Divergence 

2.14 

375 

“124 

150 


^.0 

6.57 

8.28 

.589 


Test series VII 


0.1410 

.1410 

0.0708 

.0708 

0 

0 

0 

0 

5.41 

5.41 

13-7 

13-7 

Divergence 

Flutter and divergence 

2.14 

2.14 

522 

kQl 

284 

fi 35 

292 

2kQ 


' 14.7 

4.86 

4.48 

2.64 

ei.25 

1.146 

.974 

.1410 

.0875 

0 

.08 

5.30 

13.0 

Flutter 

2.14 

486 


253 

(a) 


4.51 

1.11 

.993 

.1580 

.1217 

0 

-.26 

3.40 

11.1 

Divergence 

2.14 

427 

f 221 

195 


^3.3 

4.90 

3.26 

.766 

.1750 

.1475 

0 

-.14 

2.97 

8.1 

Divergence 

2.14 

393 

^99 

165 


*^.1 

6.16 

4.46 

.648 

.1750 

.1475 

0 

.36 

2.87 

8.1 

Flutter 

2.14 

332 

% 

118 

4.30 

2.0 

5.20 

.96 

.465 


“Not obtained. 

^Calculated flutter speedj the calculated divergence speed is 487. 
“Calculated flutter frequency. 

“struts covered with fairing simulating thin airfoil. 

*^Calculated flutter speed; the calculated divergence speed is 284 . 
®Based on calculated flutter speed. 



Vo) 

VJl 


NACA RM L55EO7 RESTRICTED 











25 - 



SfrtcHtiMd tofl steta 


(a) Airfoil- shaped body^ basic 
open tube, and supports. 


(b) Bodies of revolution used as nose 
and tail sections for the basic 
tube . 


Figure 1.- Dimensions of the various bodies and the support system. 
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(a) Open tube. 

Figure 2.- Models mounted in the Langley 4.5-foot flutter research tunnel. 



(b) Airfoil- shaped body. 
Figure 2.- Concluded, 
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Figure 4.- Plot of frequency of yawing oscillations against airspeed for 
the airfoil- shaped "body. (Table 5, test series II.) 
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DlmenslonlesB center-of-gravlty location, 

Figure 7-- Plot of speed coefficients at flutter and divergence against 
dimensionless center -of- gravity position for the open tuhe on lanfaired 
charts. (Table 5» test series VI.) 
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Symbols without flag: m = 0.l4l0 
Symbols with flag at left: m = .I 58 O 
Symbols with flag at right: m = .1750 


□ Flutter 

O Divergence 
Q Flutter and 
divergence 

A Flutter analysis 
bk Divergence analysis 


Dimensionless center-of -gravity location, 

Figure 8.- Plot of speed coefficients at flutter and divergence against 
dimensionless center-of-gravity position for the open tube on faired 
charts. (Table 5, test series VII.) 
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Dimensionless elastic-axis location, a Dimensionless dynamic pressure, q* 

(a) Speed coefficients at onset of (t) Frequency ratios at onset of 
yawing oscillations. yawing oscillations. 

Figure 9.- Speed coefficient and frequency ratio at onset of yawing 
oscillations of various bodies of revolution. (Table 3, test 
series III.) 
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